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Exponential stability of time-delay systems via new 

weighted integral inequalities 

L.V. Hien and H. Trinh 


Abstract —In this paper, new weighted integral inequalities 
(WIIs) are first derived by refining the Jensen single and double 
inequalities. It is shown that the newly derived inequalities in 
this paper encompass both the Jensen inequality and its most 
recent improvements based on Wirtinger integral inequality. 
The potential capability of the proposed WIIs is demonstrated 
through applications in exponential stability analysis for some 
classes of time-delay systems in the framework of linear matrix 
inequalities (LMIs). The effectiveness and least conservativeness 
of the derived stability conditions using WIIs are shown by 
various numerical examples. 

Index Terms —Exponential estimates, time-delay systems, 
integral-based inequalities, linear matrix inequalities. 


I. Introduction 

The problem of stability analysis and its applications to con¬ 
trol of time-delay systems is essential and of great importance 
for both theoretical and practical reasons 0 ]. This problem has 
attracted considerable attention during the last decade GI-0- 
Many important results on asymptotic stability of time-delay 
systems have been established using the Lyapunov-Krasovskii 
functional (LKF) method in the framework of linear matrix 
inequalities (LMIs) 0. It is a fact that asymptotic stability is a 
synonym of exponential stability 0, and in many applications, 
it is important to find estimates of the transient decaying rate 
of time-delay systems 0- Therefore, a great deal of efforts 
has been devoted to study exponential stability of time-delay 
systems 0-021. To derive an estimate, also referred to as 
a-stability, of the exponential convergence rate of a time- 
delay system, various approaches have been proposed in the 
literature. For example, state transformation £(£) = e at x(t) 
combined with the Lyapunov-Krasovskii functional method 
m-m, model transformation 03, constructing modified 
LKFs with exponential weighted functions 0, D3-1I91, esti¬ 
mating the Lyapunov components D3 or modified comparison 
principle 0, l20l . 

However, looking at the literature, it can be realized that 
the proposed methods in the aforementioned works usually 
introduce conservatism in exponential stability conditions not 
only on the exponential convergence rate but also on the 
maximal allowable delay and the number of matrix variables. 
Therefore, aiming at reducing conservativeness of exponential 
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stability conditions, ant important and relevant issue is to 
improve some integral-based inequalities. 

In this paper, we first propose some new weighted integral 
inequalities (WIIs) which are suitable to use in exponential 
stability analysis for time-delay systems. We show that the 
newly derived inequalities in this paper encompass both the 
Jensen inequality Gil and some of its recent improvements 
based on Wirtinger integral inequality 0, H221. We then 
employ the proposed WIIs to derive new exponential stability 
conditions for some classes of time-delay systems in the 
framework of linear matrix inequalities. Numerical examples 
are provided in this paper to show the efficiency and potential 
capability of the newly derived WIIs. 

The rest of this paper is organized as follows. In Section 
2, some preliminary results are presented. New weighted 
integral inequalities and their applications in exponential sta¬ 
bility analysis for some classes of time-delay systems are 
presented in Section 3 and Section 4, respectively. Section 5 
provides numerical examples to demonstrate the effectiveness 
of the obtained results. The paper ends with a conclusion and 
references. 


II. Preliminaries 


It can be realized in many contributions that, to derive 
the exponential estimates for time-delay systems, a widely 
used approach is the use of weighted exponential Lyapunov- 
Krasovskii functional 0. For example, a functional of the 
form 


i-O rt 

V(xt) = / x T (u)Rx{u)duds (1) 

J—r Jt-\-s 


where x is the state vector, scalars a > 0, r > 0 and matrix 
R > 0, has been used in many works in the literature D3- 
D2l The derivative of V (Xt ) is given by 


V(x t ) = tx t {t)Rx{t) 


e Q ( s ^x T (s)Rx(s)ds. 


/ t — T 


( 2 ) 


In order to generate LMIs conditions, an estimate on the 
second term of 0 is obviously needed. The problem raised 
here is how to find a tighter lower bound of a weighted integral 
of quadratic terms in the following form 


t-b 

I w (ip,a) = / e a( ' s ~ b ' ) ip T (s)Rxp{s)ds 
J a 


where a > 0 is a scalar, ip £ C([a, b], IR") and R is a 
symmetric positive definite matrix in R" XT1 , R £ S+. When 
a = 0 we write I(<p) instead of I w (ip, 0). 


Inspired from the proof of the Jensen inequality ED , we 
have the following results which referred in this paper to as 
Jensen-based weighted integral inequalities (WIIs) in single 
and double forms. 


Lemma 1. For a given matrix R £ Sf, scalars b > a, a > 0, 
and a function ip £ C([a, b], R"), the following inequalities 
hold 

r b , T . rb 


(/ <p(s)dsj R^j p(s)ds (3) 


0 a(u—b) T 


ip ( u)Rip(u)duds 


qA . rb rb rO rO 

— \ J J F{ u )dudsj Ry J J ip(u)dudsj , 


> 


(4) 


where 7 fc = e‘ 


— pOt(b—a) _ 


E k 

.7=0 


fb-ay 


, k > 0. 


J j=0 o'- 

Proof: By taking integral of the inequality 

> 0 we obtain 


e “( s — fe )v5 T (s)i?(p(s) ip T {s) 


ip(s) 


gOt(b—s) R —l 


I w {ip,a) f b ip T (s)ds 
Ja P{s)ds p(a)R~ 1 


> 0 


which implies © by Schur complement. The proof of © is 
similar and thus it is omitted here. ■ 

Remark 1. Obviously > e -“( h - a ) for all a > 0, 6 > a. 
Therefore, © gives a new lower bound in comparison to the 
common estimate I w (ip,a) > e -“( fc -“) Ifp). 


Remark 2. When a approaches zero the previous inequalities 
lead to the Jensen inequality in single and double form, respec¬ 
tively. More precisely, from the fact that lim Q _> 0 + 
we readily obtain the following results 


j ip T (s)Rip(s)ds > ^ J ip(s)ds^ r(^ J ip(s)d. 

p; 

rb rb 

/ / ip T (u) Rip(u)duds 

J a J s 

(/ J tp(u)duds^J r(^ J J ip(u)duds' S J. 


(b-a) 2 


III. New weighted integral inequalities 


In this section, some new weighted integral inequalities are 
derived by refining ©, ©. In the following, let us denote 


= I w (v,u) - ~{j v( s ) ds ) R(j V(s)ds ) 

as the gap of ©. By refining © we find a new lower 
bound for J w (ip,a) other than zero. First, let us introduce 
the following notations for given scalars b > a, a > 0, and 

¥>£C7([a,&],R n ) 


£ = b — a, A a 


To 


(l+7o)l 2 

To 


Li=[l -^],C = col|y ip(s)ds, j y <p(u)duds 


By using the Taylor series expansion of exponential function, 
it can be verified that A a > 0 for any a > 0. We also 
use the notion of Kronecker product A <g) B for matrices 
A £ IR nxrra ,JJ £ R ryXr . For more details about the Kronecker 
product, we refer the readers to 11231 . 


Lemma 2. For a given nxn matrix R > 0, scalar a > 0, and 
a function ip £ C([a, 6], R n ), the following inequality holds 

JS,(<P,a) > -CiL^L^RX (7) 

Po 

where p Q = A a = 3) [ 7 g - (af) 2 e at ]. 


Proof: For any ip £ C([a , b ], R"), we define an approxi¬ 
mation function ip £ C([a, 6], R") as follows 

ae a(b-t) rb 

ip(t) = p{t) -/ ip(s)ds + h(t)x (8) 

To Ja 

where h(t) is a real valued function on [a, b] and \ G lf& n is 
a constant vector which will be defined later. For brevity we 
let w(t) = e a ( b-t ) and predefine h(t) = w(t)p(t), where p(t) 
belongs to P&, the set of polynomials of order less than k. A 
prior computation gives 


Jw(y> a ) = Jw(T,Oi) + J w {h)x T R\ 

rb rb 

— 2 — / h(s)dsx T R / (f(s)ds 
70 J a J a 

( rb rb rb 

p(a) / ip(s)ds + p'{a) / / ip(u)duds 

J a J a J s 

rb rb rb \ 

+ / p"( s ) / / <p{v)dvduds J, 

J a J s J u J 

where J w (h) = J b w~ 1 {s)h 2 (s)ds - J* h(s)ds ) 

Now, for any p £ Pi which we can write p{t) = Co + cit, 
ci f 0. Then 


(9) 


U, i , p{d)e at — p(b) , Ci 7o 

h(s)ds =-1- 

a or 

i 2 (s)ds = eQ V( a ) -P 2 (.b) 
a 

2ci (e ae p(a) - p(b)) 2cfy 0 
a 2 a 3 

and thus J w [h) = 4r c i- From © we obtain 

a ajo 

( 10 ) 

By Lemma [I] J^,{ip, a) > 0 which leads to 

> -^^x T Rx+‘^-^X T R(L 1 ®I n )(. (11) 

a a 7 o 

Hereafter, we will denote by 7 Z(J‘j J (ip, a )) the right-hand side 
of (11 111 . Now we define vector x i n the form X = 
where A is a scalar, then 

7 im<p,a)) = — ^2—A - A a A 2 )c r (^i ® 7?)C 

a \ 70 / 
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A. Systems with discrete and distributed constant delays 


The function f( A) = 2 —A — A a X 2 attains its maximum , 11 .. 
at A = A 7l ^ o . Then it follows from (jTTJ that J^(ip,a) > 
( T (LjL i 0 i?)£ which completes the proof. ■ 

Remark 3. It is interesting that estimate 0 does not depend 
on the selection of first-order polynomial p G Pi- In other 
words, inequality 0 can be derived from 0 for any function 
h(t) of the form (co + ci ^ 0. Of course, when 

Ci = 0 then 0 leads to ®. 

Remark 4 . By repeating the proof of Lemma [2] with the 
approximation 

a 2 w(t) f b I b 

fi(t) = tp(t) -/ / p(u)duds + w{t)p{t)x (12) 

7l J a Js 

where w(t ) = and p G Pi then ® leads to double 

WII formulated in the following lemma. 


Lemma 3. For a given nxn matrix R > 0, scalar a > 0, and 
a function ip G C([a, 6], P"), the following inequality holds 

e a ^ u ~ b ^ p T (u)Rip{u)d U ds (13) 

> — C T (L^L 0 ®f?)C+— ( t (LTL 2 ®R)( (14) 
7i Pi 



where £ = col |/ a 6 (p(u)duds, f b J b (p(v)dvduds'^, 


L 0 = [1 0], L 2 = [1 -^] 

TD 271 _ \a£+(at— l)7o]l 2 

O' q,2 7i 


Pi = 


I Q!7i 

V 72 


)■ 


B n and 


Remark 5. The following facts can be found by using Taylor 
series of the exponential function 


lim — 

a— >0+ po 

li m * 

a—>0+ pi 


3 

b — a ’ 


lim L\ = Li 
a— >o+ 


[1 


—1 

b-oj i 


4 

{b -a) 2 ’ 


lim L 2 = L 2 

a—>-0+ 


[i 


-3 
b—a. 


Therefore, when a approaches zero we obtain the following 
results which are the same as those derived by the Wirtinger 
inequality in single and double form EH, EH! 

[ p T (s)Rtp(s)ds > — C t {LqL 0 + 3L^Li) <g>i?£ (15) 

Jn 0 CL 


ip T (u)R(p(u)duds 

> n 2 , 2 C T {LqL 0 + &L 2 L 2 ) 0 R(. (16) 

(6 — a) z 

Remark 6. The results obtained in Lemma \2\ and Lemma [5] 
can be extended using 0 and (fl2l> where h(t) = e a( - b *1 p(t) 
and p(t) belongs to the set of higher order polynomials, and 
then new lower bounds for 0. ® can be derived. This will 
be addressed in future works. 



IV. Exponential estimates for time-delay systems 

This section aims to demonstrate the effectiveness of the 
newly weighted integral inequalities proposed in this paper 
through applications to exponential stability analysis for two 
classes of time-delay systems. 


Consider the following time-delay system 

x(t) = A 0 x(t) + Ayxit — h) + A 2 f x(s)ds, t > 0, 

J t—h 

x(t.) = fit), t g [-M], 

(17) 

where Aq,Ai,A 2 G P" x " are given constant matrices, 
h > 0 is known time-delay, f G C([—h, 0], IP") is the initial 
condition. 

We recall here that q. m, for a given a > 0, system 
(ITTI) is exponentially stable with convergence rate a if there 
exists /? > 0 such that any solution x{t, <j>) of lITTl) satisfies 
\\x(t,<t>)\\ < /3||0l|e _crt , Vf > 0. 

Let a G P" x4 " defined by a = [0„ x (j_i)„ I n 0 rax ( 4 _ i ) n ], 
i = 1,..., 4. We denote A = Aoe i + Aie 2 + A 2 e 3 and the 
following matrices 



ei 


A 

F 0 = 

e 3 

, Fi = 

ei - e 2 


f4_ 


he 1 - e 3 _ 


F 2 


ei - e 2 
he i - e 3 


F » = 


he i - e 3 
h 2 / 2ei — e 4 


n 0 = if PF 1 + FT PF 0 + a if PF 0 , 

IIi = eiQe i — e~ ah e 2 Qe 2 + A T (hR + h 2 /2Z)A, 

H 2 = —(e 4 — e 2 ) T R{ei — e 2 ) 4-if (if Li 0 R)F 2 , 

70 Po 

2 a 2 

n 3 = —(hei - ef) T Z(he\ - e 3 ) + -^-if (L 2 L 2 0 Z)F 3 

71 Pi 


where L\,L 2 and 70 , 71 , po, pi are defined in ©, ©, 0. (Q® 
with a = —h, 6 = 0. 

Theorem 1. Assume that, for a given a > 0, there exist 
symmetric positive definite matrices P G |R 3,l><3rl ; Q,R,Z G 
P" x " satisfying the following LM1 


n 0 + iii ~ 1I2 — n 3 < 0. (18) 


Then system o is exponentially stable with a convergence 
rate a = a/ 2. 

Proof: Consider the following Lyapunov-Krasovskii func¬ 
tional 


V(x t ) = x T (t)Px(t) + f e“^ s t?> x T (s)Qx(s)ds 

J t—h 

/•0 pt 

+ / e a (“ _t ) x T (u)Rx(u)duds (19) 

J—h J t-\-s 

n t 

e a< ' e ~ t )x T (d)Zx(9)duds 

-\-U 


where x(t) = 


(t) = col |x(f), f*_ h x(s)ds, f*_ h fl x(u)dud s y 


It follows from ( IT9l that V(x t ) > A m i n (P)||a;(t)|| 2 . Taking 


3 


















derivative of V(xt) along trajectories of (ITTI) we obtain 

V(x t )+aV(x t ) = x T (i)( n o+ IIi)x(f) 

(20 ) 

( e a ( u_t ) ± T (u)Zx(u)duds 

t-\-s 

where 

X(f) = col {x(t),x{t - h),f t t _ h x(s)ds,f° h f t t +s x(u)duds}. 

By applying Lemma U and Lemma [3] to the first and the 
second terms in (l 20 l >. respectively, we then obtain 

v(x t ) a aV(x t ) < x T (t)( n 0 + n 4 - n 2 - n 3 ) x (i). (2i) 

It follows from (IT 8 ] > and ( |2T1 > that V(xt) + aV(xt) < 0 , 
which yields V(xf) < V((j>)e~ at . This leads to ||x(f, 0)|| < 

\j' \ V ^(P ) e ~ a ^ 2t • The P ro °f i s completed. ■ 

Remark 7. Note that the exponential transformation z(t) = 
e at x(t),cr > 0, in general, is not applicable to access the 
exponential stability of system (Ell because it leads to a 
system with time-varying coefficients. When 3l 2 = 0, using 
the aforementioned transformation, system O becomes 

z(t) = (A 0 + al n )z(t) A e ah A 1 z(t - h). (22) 



In many works found in the literature, in order to get expo¬ 
nential estimates for system (O (with A 2 = 0 ), it was trans¬ 
formed to ( 122 b first and then asymptotic stability conditions for 
(122b were proposed. However, this approach usually produces 
conservatism in exponential stability conditions due to the fact 
that the exponential stability of (O (with A 2 = 0) is just 
equivalent to the boundedness of <0 which is less restrictive 
then asymptotic stability. Differ from those, and as discussed in 
Section 2 in this paper, we here propose an improved approach 
used in exponential stability analysis for time-delay systems by 
employing our new weighted inequalities derived in Lemma \2} 
and Lemma [i] in this paper. 


B. Systems with interval time-varying delay 


Consider a class of linear systems with interval time-varying 
delay of the form 

( x(t) = Ax(t) A A d x(t - h(t )), t > 0 
\x(t) = <f>(t ), t G [-h 2l 0] 


T 0 = col{.4, ei - e 2 , e 2 - e 4 }, 

Ti = coljei - e 2 ,h\(e\ - e 5 )}, 

T 2 = col{e 2 - e 3 , e 2 + e 3 — 2e 6 }, 

T 3 = col{e 3 - e 4 , e 3 + e 4 — 2e 7 }, A = [T 2 TJ] t , 
fi 0 (/i) = T(h) r PT 0 + T x PT(h) + aT(h) T PT(h), 
fit = ejQ 4 ei - e~ ahl e2Qie 2 

+ e~ ahl e^Q2e2 - e~ ah2 ejQ 2 ei, 
n 2 = A T (h\Ri A h 2 12 e ahl R 2 ) A, 

^3 = —r—(ei - e 2 ) T f?i(ei - e 2 ) + <g> Ri)Y u 

7o Po 

7o = e ahl — 1, 7i = 7o — ahi, L 1 = [1 - 0 : 70 / 71 ], 

Po = %, [70 - {ahi) 2 e ahl ] , hi 2 =h 2 -h\. 

7i 

Theorem 2. Assume that there exist symmetric positive def¬ 
inite matrices P G j^ 3 nx 3 n^ q. p>. g R" xra , i = 1,2, and 
a matrix X G R 2 "x 2n such that the following LMIs hold for 
h G {hi, )i 2 } 




(24) 


Q(h) = n 0 (h) A fli + n 2 - Vs - e"“ ftl 2 A T nA < 0 , (25) 


where f ? 2 = diag{i? 2 , 3f? 2 }. Then system (123b is exponentially 
stable with a convergence rate a/2. 


First, we need the following lemmas. 

Lemma 4. If fl(hi) < 0 and f l(h 2 ) < 0 then Ll(h) < 0, 
Wi G [hi, h 2 \ 

Proof: It is obvious that = 2aT T PT > 0, 

where T = [0 7raX 27! (e 6 —e 7 ) T ] T . Therefore, f 1(h) is a convex 
quadratic function with respect to h. This completes the proof. 


Lemma 5 (Improved reciprocally convex combination j24l ). 
For given symmetric positive definite matrices R\ G IR" X? \ 
li> G IR mxm , if there exists a matrix X G IR nxm such that 

1 ^ > 0 then the inequality 


i ?2 


r 1 




0 


0 t^R2 

holds for all 6 G (0,1). 


> 

R i X 


* R'2 


where A, A d G IR" xn are given constant matrices, h(t) is 
time-varying delay satisfying 0 < hi < h(t ) < h 2 , where 
hi, h 2 are known constants involving the upper and the lower 
bounds of time-varying delay. 

Let ei [0nx(z— l)n In ^nx(7—i)n]-> t 1,2, ...,7. We 
denote A = Ae i + A d ez and 


Xi(f) = col ■ 


f 

x(t) 



x(t — hi) 



x(t — h(t)) 

1 


x(t — ft 2 ) 

- 


7u /t-hi x ( s )ds 
h(t)-h i ft.-h(t) x ( s ) ds 

h^W) St-™ X ^ ds . 


T (h) = col{e 4 , hie 5 , (h - hi)e 6 + (h 2 - h)e 7 }, 


Proof: Inspired from ED, we now consider the following 


LKF 


V{x t ) = Xo (t) p Xo(t) A [ e“ (s t) x T (s)Q 1 x{s)ds 
J t—hi 


rt—h i 


+ 


o a ( s - t )r T 


x 1 ( s)Q 2 x(s)ds 


(26) 


J t — h 2 

/■o H 

Ahi / £«(«-*) + T 

J — h\ Jt-\~s 

p — h\ pt 


x T (u)R\x(u) duds 


Ah 


12 


= a(u+7ii-t) • T 


h -2 Jt + s 


x 1 ( u)R 2 x(u)duds 


4 



















where xo(t) = col{z(i), J*_ hi x(s)ds, f*_ £ x(s)ds}. 

It follows from (l26t that V(x t ) > A m i n (P)||;r(f)|| 2 . In re¬ 
gard to the fact xo (t) = G 0 (h)xi(t) and j|xo(f) = GiXi(t), 
the derivative of (126} along trajectories of ( 1231 ) gives 

V(x t ) + aV(x t ) = Xi(t) (do(h) + fil + n 2 )xi (t) 

- hi [ e a( ' s ~ t) x T {s)Rix(s)ds 

Jt-h i y Lt > 

-h 12 f 1 e a{s+hl - t) x T {s)R 2 x(s)ds. 

J t—h 2 

By Lemma [2] we have 

-hi f e a( ' s ~ t ' , x T (s)R 1 x(s)ds <-Xi(t)^3Xi(t)- ( 28 ) 

J t—h\ 


Corollary 1. System (123b is asymptotically stable for any 
delay h(t ) E [ft.i, / 12 ] if there exist symmetric positive definite 
matrices P E R 3 " x3ri , Q t . R t g |R" X ", * = 1,2, and a matrix 
X E fj 2 nx 2 n sa ti s fying (124b and the following LMI holds for 
h G {hi, h 2 } 

n o (h) + <f>i~$ 2 -A T IlA<0, (31) 

where 

n 0 (h) = T(h) T PT 0 + T^PT(h), 

$1 = e-iQiei - e 2 Qie 2 + e 2 Q 2 e 2 - ej<3 2 e4 

4- (/i 2 i?i + h\ 2 Ri) A, 

<E> 2 = Tjdiag{P 1 ,3i? 1 }T 4 , 

T4 = col{ei — 22, ei + e2 — 2es}. 


Next, by splitting 

pt—hi 


x(s+h!-t) - T 


't—h 2 


ar (s)i?2^(5)ds 


/»t—/ii 


J t—h(t) 

pt—h(t ) 

/ 

' t—h 2 


gocls+tn- t )x T {s)R 2 x{ S )d S 


p a(s+hi-t) ■ T 


x 1 ( s)R 2 x{s)ds 


the second integral term of (127} can be bounded by (ITT} and 
Lemma [4] as follows 

—hi 2 f e a( ' s+hl ~ t ^x T {s)R 2 x(s)ds 

J t—h2 

hi 2 e~ ahl2 T . ,, hi 2 e~ ahl2 - - 


< -■ 


h(t) — h 


-X2 {t)R2X2(t) - 


< —e 


ahi2 

X2 (t) 

T 


X3 (t\ 


ahi2 

X2 (t) 

T r 


X3 (t) 



n 12 E> 
h(t)—hi 2 

0 

#2 A" 

* R 2 


h 2 - h(t) 
0 


X'3 (t)R2X3(t) 


n 12 D 
h2~h(t) 2 

X2(f) 

X3 (i) 


X2(f) 

X3(i) 


where x 2 (f) = T 2 Xi(i) and X3<7) = Y 3 xi(t). Note that the 
previous inequality is still valid when h(t) tends to hi or h 2 . 
This leads to 

nt — hi 


-hi 


0 a{s+h a -t) -T 


t t—h 2 


x 1 ( s)R 2 x(s)ds 
< -e~ ahl2 Xi (t)A T HA X i(t). 


(29) 


Combining (|27} -(|29} we then obtain 

V(x t ) + aV{x t ) <xl(t)n(h)xi{t). (30) 

By Lemma[4] ([25} implies that 12 (h) < 0 for all h G [hi, h 2 \. 
Therefore, if (125} holds for h = hi and h = h 2 then, from 
(l30b . V(xt) + aV(xf) < 0 which concludes the exponential 
stability of (123} with guaranteed decay rate er = a/2. The 
proof is completed. ■ 

Remark 8. When a approaches zero, by Remark [5] and 
Theorem\2\ we obtain the same asymptotic stability conditions 
for system (ED) derived from improved Wirtinger’s inequality 

SZfl. 


V. Examples 


Example 1. Consider the following system ©, ll22l 


x{t) 


0.2 

0.2 


0 

0.1 


x{t) + 


-1 

-1 


0 

-1 



(32) 


An eigenvalue analysis a shows that m remains stable 
for a constant delay in the range h G [0.2,2.04]. By the 
Wirtinger-based inequality approach. Theorem 6 in © and 
Theorem 1 in m guarantee the asymptotic stability of (132} 
for h in the interval [0.2,1.877] and [0.2,1.9504], respectively. 
In Theorem Q] we fix a at 0.0002, it is found that (ITS} is 
feasible for h G [0.2,1.9778]. This clearly shows a reduction 
of conservatism of Theorem Q] 


Example 2. Consider system ([17} with the matrices taken 
from the literature 


A q = 


0 

-2 


1 

0.1 ’ 


A i 


0 0 
1 0 ’ 


A 2 


0 0 
0 0 


It is surprising that, for this system, the exponential stability 
criteria proposed in GH2), especially 03, are not feasible 
for any h > 0, a > 0. By using the notations given in 

|4l , we have A4 = A/f + |Ai| = q n \ ■ Thus, for any 


= (fi,^2) T G 


3 0.1 

^ 2 , vi > 0,^2 > 0, Mv = 


> 0. 


v 2 

3iq + v 2 \ 

This shows that the stability conditions proposed in |4j and 
EO) are not applicable to this system. Note also that, since 
Re(eig(Ao + Ai)) = 0.05 > 0, the delay-free system 
is unstable and the results to access exponential stability 
of this system are much more scarce. By employing the 
Wirtinger-based integral inequality, a significant improvement 
of the asymptotic stability conditions was provided in ©. To 
compare with our approach, we apply Theorem 6 in © to 
(l22b . The obtained results in Table 1 show that, thanks to 
our new weighted integral inequalities proposed in Lemma [2] 
and Lemma [3] Theorem Q] in this paper gives significantly 
better results. The simulation result presented in Figure 1 is 
taken with h = 1.6, er = 0.045, incremental step s = 10~ 4 
and initial condition </>(t) = [2 — l] r . It can be seen that 

the state trajectory z(t) = e at x(t) is bounded, and thus. 
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TABLE I 

Decay rate a for various h in Example 5.2 


h 0.3 

0.5 

0.8 

1.0 

1.5 

1.6 

NoDv 

~m (L 097 T 

0.1905 

0.2936 

0.2766 

0.0175 

- 

3 ri z + 2 n 

Theorem □ 0.0971 

0.2095 

0.4195 

0.4978 

0.1039 

0.045 

6 n 2 + 3n 


NoDv: Number of Decision variable 



Fig. 1. Trajectory z(t) = e 0 045t a:(t) with h = 1.6. 

x(t) exponentially converges to the origin with decay rate 
cr = 0.045. 


the exponential stability criteria proposed in a. Ea based on 
positive system approach also cannot access the exponential 
stability of the system. In Ifl8l . some integral equalities were 
used to overcome the conservative estimates. However, when 
manipulating the derivative of the Lyapunov-Krasovskii func¬ 
tional, all the integral terms were abandoned (see, Eq. (10) 
in [HI). This leads to the fact that the proposed conditions in 
[|T8l are very conservative. We apply the main theorem in lfl8l 
to this example, the obtained results for hi = 1 and various 
h -2 are listed in Table 3. In m, exponential convergence rate 
of solutions was derived by estimating the maximal Lyapunov 
exponents. By Theorem 3 in fl4l . the exponential convergence 
rate a £ (0, A*], where A* is the unique positive solution 
of equation A + 0.0087e A,t2 = 0.2707. We apply Theorem 
U in this paper for hi = 1 and various /#2 • The obtained 
results are exposed in Table 3. Clearly a significant reduction 
of conservatism is delivered by Theorem [2] This shows the 
effectiveness of our approach. The simulation result presented 
in Figure 2 is taken with h(t) = 1 + 5|sin(f)|, a = 0.2546 
and incremental step s = 10” 4 which illustrates our theoretical 
results. 


Example 3. Consider an active quarter-car suspension system 
with control delay introduced in l25l 

x(t) = Ax{t) + Bu(t-h(t)), 
y(t) = Cx(t), 

where x(t) £ [R 4 is the state, u(t) is the control input, y(t) is 
the output and 



' 0 

0 

1 

A = 

0 

k.. 

0 

0 

0 

c s 


771 s 

k s 
- m u 

rri £ 

kt _£s_ 

m u m u 


" 0 ' 


'1' 

B = 

0 

1 

, c = 

1 

1 

0 


m u . 



The following parameters are taken from |[25l . 



Fig. 2. Trajectory e 0 2546t x(t) with h(t) = 1 + 5| sin(i) |. 


TABLE II 

Quarter-car model parameters 


Example 4. Consider system (l23l) with 


m s 

m u 

ks 

kt 

C s Ct 


0 

1 

, A d = 

0 

0.1 

973kg 

114kg 

42720N/m 

101115N/m 

1095Ns/m 14.6Ns/m 

A — 

-10 

-1 

0.1 

0.2 


A static output feedback controller is proposed as u(t) = 
Ky(t), where K is the controller gain. The closed-loop system 
is then in the form of (l23l > with A,i = BKC. For illustrative 
purpose we consider K = 1. It is noted first that, in this case. 


This example has been taken from l26l . The obtained results 
by Corollary Q] as listed in Table 4. These results again show 
the effectiveness of our approach in proposed in this paper. 
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TABLE III 

Decay rate <t for hi = 1 and various h 2 


h2 

2 

3 

4 

5 

6 

NoDv 

[18| 

0.1338 

0.1316 

0.1300 

0.1290 

0.1267 

50.5n 2 + 6.5n 

□3 

0.2562 

0.2522 

0.2473 

0.2416 

0.2351 

- 

Theorem [2] 

0.2690 

0.2672 

0.2644 

0.2603 

0.2546 

10.5n 2 + 3.5n 


TABLE IV 

Upper bound of h 2 for various hi = 1 


hi 

0 

0.3 

0.7 

1 

2 

NoDv 

"HU 

0.55 

0.77 

1.17 

1.47 

2.47 

50.5n 2 + 6.5n 

03 

1.35 

1.64 

2.02 

2.31 

3.31 

9.5n 2 + 5.5n 

(3D 

1.64 

2.13 

2.70 

2.96 

3.63 

21 n 2 + 6 n 

Corollary 

1.88 

2.18 

2.53 

2.81 

3.78 

10.5n 2 + 3.5n 


VI. Conclusion 

In this paper, new weighted integral inequalities (WIIs) have 
been proposed. By employing WIIs, new exponential stability 
criteria have been derived for some classes of time-delay sys¬ 
tems in the framework of linear matrix inequalities. Numerous 
examples have been provided to show the potential of WIIs 
and a large improvement on the exponential convergence rate 
over the existing methods. 
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